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ABSTRACT 

Baryon acoustic oscillations (BAOs) and weak lensing (WL) are complementary probes of cosmology. 
We explore the distance and growth factor measurements from photometric BAO and WL techniques 
and investigate the roles of the distance and growth factor in constraining dark energy. We find for 
WL that the growth factor has a great impact on dark energy constraints but is much less powerful 
than the distance. Dark energy constraints from WL are concentrated in considerably fewer distance 
eigenmodes than those from BAO, with the largest contributions from modes that are sensitive to 
the absolute distance. Both techniques have some well determined distance eigenmodes that are not 
very sensitive to the dark energy equation of state parameters wq and w a , suggesting that they can 
accommodate additional parameters for dark energy and for the control of systematic uncertainties. 
A joint analysis of BAO and WL is far more powerful than either technique alone, and the resulting 
constraints on the distance and growth factor will be useful for distinguishing dark energy and modified 
gravity models. The Large Synoptic Survey Telescope (LSST) will yield both WL and angular BAO 
over a sample of several billion galaxies. Joint LSST BAO and WL can yield 0.5% level precision on 
ten comoving distances evenly spaced in log(l + z) between redshift 0.3 and 3 with cosmic microwave 
background priors from Planck. In addition, since the angular diameter distance, which directly affects 
the observables, is linked to the comoving distance solely by the curvature radius in the Friedmann- 
Robertson- Walker metric solution, LSST can achieve a pure metric constraint of 0.017 on the mean 
curvature parameter Oj, of the universe simultaneously with the constraints on the comoving distances. 
Subject headings: cosmological parameters — distance scale — gravitational lensing — large-scale 
structure of universe 



1. INTRODUCTION 

The luminosity distance-redshift relation measured 
from type la supernovae (SNe) has provided a vital piece 
of evidence for an accelerated cosmi c expansion in a 
Friedman- Robert son- W alker universe (|Riess et al.lll99ct 
iPerlmutter et a l. 1999). Such an acceleration suggests 
that either the universe is dominated by a smooth dark 
energy component or we need an alternativ e framework 
of gravity to accoun t for the SN data (e.g., iDvali et al.l 
2000; Def faveti l200l . Although different answers have 
vastly different theoretical implications, they can only 
be distinguished indirectly, for example, through mea- 
surements of the cosmic distance scale as a function of 
redshift and the amplitude of the density fluctuations as 
a function of redshift, i.e., the growth factor-redshift re- 
lation. One may further derive a phenomcnological dark 
energy equation of state (EOS) as a function of redshift, 
w(z) 1 from the distances and growth factors and contrast 
it with predictions from theory. 

Besides SNe, baryon aco ustic oscillations (BAOs) in 
the g a laxy power sp e ctrum (Eiscuste in. Hu. fc Tegmarkl 
1998t ICoorav et all 120011 iBlake fc Glazebrookl 120031: 



Hu fe Haiman 



Linder 



ing (WL, IHu fc Tegmarl< 



a nd wea k lens- 
, iMeniej [19991: 
Bartelmann fc Schneider! 120011: iHutererl 120021 : 
Refregierl 120031: ISong fc Knoxl |2004|) can also mea- 
sure the distance (ISeo fc Eisenstein 20031: ISone 120051 : 



iKnox. Song, fc TvsoiJl2006at iZhan fc Knoxll2006ah . The 
sound horizon at the last scat tering surface determines 
the BAO scale (iPeebles fc Yulll970l: iBond fc Efstathioul 
[1981 lHoltzmanlil989t ). which is fixed m comoving space, 



so that it can be used as a standard ruler to measure the 
comoving angular diameter distance Da(z) and, with 
spectroscopy, Hubble parameter H{z). With precise 
calibration of the sound horizon from cosmic microwave 
background (CMB) observations, the BAO technique 
is able to measure the absolute distance, whereas the 
SN technique measures the relative distance because 
of the degeneracy between the SN intrinsic luminosity 
and the Hubble constant Hq. The WL technique does 
not depend on power spectrum features; its ability to 
measure distance originates from its geometric lensing 
kernel and the sensitivity of the shear p ower spectrum 
amplit ude to the absolute distance (e.g., I Zhan fc Knoxl 
2006^ hereafter ZK06b). 

The WL shear signal is a direct indicator of the gravita- 
tional potential of all matter, luminous or not, so the WL 
technique has the advantage of being able to measure the 
growth factor of the large-scale structure directly. For 
the BAO technique 1 , the growth factor is largely degen- 
erate with the galaxy clustering bias, though it could be 
determined from the shape of the galaxy power spectrum 
if one knew precisely how the amplitude of the linear 
matter power spectrum affects the slope of the nonlinear 
matter power spectrum and how the galaxy bias behaves 
in the mildly nonlinear regime (see Section l3.2p . 

The relative importance of the distance and growth 
factor in constraining the da rk energy EOS with 
WL has been discussed by lAbazaiian fc Dodelsonl 

1 Although the BAO features are not related to the growth fac- 
tor, the BAO technique, being a power spectrum or correlation 
function analysis, can potentially measure the growth factor. 
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(l200l. iSimpson fe Bridle! (I2005I . hereafter SB05), 
Zhang. Hui fc Stebbinsl (1200,1 hereafter ZHS05), and 



Knox et all (|2006al) . but their results are apparently 
inconsistent: SB05 and ZHS05 argue that the dis- 
tance and growth fact o r are equal ly powerful whereas 
lAbazajian fc Dodelsonl (|2003l ) and iKnox et all (|2006ah 
show that the growth factor is much less powerful than 
the distance. This inconsistency is found to be caused 
by the parameter space of each analysis: with no other 
parameters, the distance and growth factor place com- 
parable constraints on the dark energy EOS parameters, 
but the growth factor becomes less powerful when more 
parameters are allowed to float (ZK06b). Because the 
WL measurements of the distance and growth factor are 
entangled, ambiguity arises on the exact meaning of "us- 
ing only the distance (or growth) information" without a 
reconstruction of distances and growth factors (actually, 
their errors). In this paper, we give a detailed account 
of estimating the distance and growth errors and their 
roles in constraining dark energy. 

Throughout this paper, we assume a low-density cold 
dark matter (CDM) universe with the following pa- 
rameters: the dark energy equation-of-state parameters 
wo and w a as defined by w (z) — wq + w a z(l + z)^ 1 
(Chcvallier & Polarski 2001), the matter density u> m , the 
baryon density u>b, the angular size of the sound horizon 
at the last scattering surface 9 S , the equivalent matter 
fraction of curvature flk, the optical depth to scattering 
by electrons in the reionized inter-galactic medium, r, the 
primordial helium mass fraction Y p , the spectral index 
n s and running a s of the primordial scalar perturbation 
power spectrum, and the normalization of the primor- 
dial curvature power spectrum at k = 0.05 Mpc -1 . 
For reasons stated in Section 13. 1\ only a subset of these 
parameters are retained when distance and growth pa- 
rameters are present. T he fiducial model is taken from 
the 3-year WMAP data (jSpergel et alj 120071) : (w , w a , 



UJ m , UJ b , 9 S , fit, T, Y, 



pi 



««> A 2 R) 



-1, 0, 0.127, 0.0223, 



0.596°, 0, 0.09, 0.24, 0.951, 0, 2.0 x 10~ 9 ). The reduced 
Hubble constant h — 0.73 and the present equivalent 
matter fraction of dark energy fix = 0.76 are implicit 
in this parametrization, meaning that either one of them 
can replace 9 S or any parameter that affects 9 S . 

The rest of the paper is organized as follows. Section [2] 
provides a brief introduction to the power spectrum 
analysis and error forecast with BAO and WL. For the 
investigation that follows, we assume a photometric 
redshift (photo-z) survey base d on the Large S ynop- 
tic Survey Telescope 2 (LSST, llvezic et all 120081 ) and 
consider only photo- z BAO and WL tomography. We 
compare BAO and WL constraints on the distance 
and growth factor in Section [3] and examine the roles 
of distance and growth for studying dark energy in 
Section |U The conclusions are drawn in Section [5] Since 
our main goal is to understand the intrinsic properties 
of the distance and growth factor measured from the 
two techniques and their applications in dark energy 
studies, we neglect systematics in Sections [3] and |4j 
To be more realistic, we include in Section [5] a brief 
discussion of the constraints on the distance and growth 
factor from a joint analysis of LSST BAO and WL 
with conservative estimates of systematic uncertainties 

2 See http://www.lsst.org 



in photo-zs and power spectra. For more detailed 
accou nts of var io us sys te matics, see iHuterer fe Takadal 



Huterer et al.l (I2006D: 
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2. WEAK LENSING AND BARYON ACOUSTIC 
OSCILLATIONS 

The projected lensing potential <j){0) transforms a 
small displacement A/3 in the source plan e into a small 
displacement A6 in the imag e plane via ([Kaiser 1992; 
iBartelmann fe Schneidei1l200lh 



-72 
K + 71 



AO, 



where the convergence k = V 2 0/2, the shear components 



7l = (d 2 (f>/d9 2 - d 2 (f>/d9 2 )/2 and 72 = d> 



89 2 . The 



shear 7 = 71 + 172 represents the distortion of an image 
without magnification, so that it can be inferred from 
the average ellipticity of galaxies within an appropriate 
window under the assumption that galaxies are randomly 
oriented in the absence of lensing. One may decompose a 
shear map into E-modes and B-modes, but, to first order, 
gravitational lensing from density fluctuations does not 
induce B-modes. Therefore, we only consider the E-mode 
shear statistics here. 

The angular power spectra of the shea r 7(6?) and galaxy 
number den sity n{6) can be written as (|Hu fe Jainll2004t 
Zhan 2006) 



Pi 



XY 



dz H{z)D A (z)W l x (z)Wj {z)Aj(k; z), 

(1) 

where lower case subscripts correspond to the tomo- 
graphic bins, upper case superscript labels the observ- 
ables, e.g., X = g for galaxies or 7 for shear, A 2 (k;z) 
is the dimensionless power spectrum of the density field, 
Da (z) is the comoving angular diameter distance of red- 
shift z as viewed from redshift [the same as the co- 
moving distance, D(z), between redshift and z in a flat 
universe], and k — £/Da(z). In the linear regime, the 
power spectrum is scaled by the growth factor G(z) 

A 2 (fc;z)=G 2 (z)A 2 (fc;0)/G 2 (0). 
The window functions are 



W?{z) = b{z)- 



W?(z)- 



3 tt m H 2 Da(z) 
2 H(z) ac 



dz 



, rii( z') D A (z,z') 
i Da(z') ' 



where b(z) is the linear galaxy clustering bias, Q m — 
Um/h 2 , and Da(z,z') is the comoving angular diameter 
distance of redshift z' as viewed from z [the same as the 
comoving distance, D(z, z'), between redshift z and z' in 
a flat universe]. The galaxy redshift distribution rii(z) 
in the ith tomographic bin is an average of the under- 
lying three-dimensional galaxy distribution over angles, 
and the mean surface density fii is the total number of 
galaxies per steradian in bin i. The binning for WL need 
not be the same as that for BAO. The irreducible occur- 
rence of the Hubble parameter H(z) in equation ([1]) for 
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BAO is due to the fact that the galaxy number density 
rii(z) is measured in redshift not in distance. 

Equation JT]) is applicable to the shear power spectrum 
for WL (X = Y = 7), galaxy power spectrum for BAO 
(X = Y = g), and galaxy-shear power spectrum (X = g 
and Y = 7). Although we focus on gaining a better un- 
derstanding of the BAO and WL techniques by exploring 
their differences, we must emphasize that a joint analysis 
of BAO and WL with all the t hree types of power spectra 
(|Hu fe Jainll200l [Zhar][2006l) is far more powerful than 
either technique alone. 

We para metrize the underlyin g galaxy redshift distri- 
bution as (|Wittman et al.ll2000f ) 

n(z) oc z a exp [-[z/z*f] 

and adopt the values a = 2, z* = 0.5, and j3 = 1 with a 
projected galaxy number density of n to t = 50 per square 
arcminute for LSST. The galaxy di stribution rij (z) in the 
ith b in is sampled from n(z) by (|Ma et al.l l2006: Zhan 
[2f)0fih 

rii(z) = n(z)V(z^ ti ,Zp/,z), 

where the subscript p denotes photo- z space, z® i and z® i 
define the extent of bin i, and V(a, 6; z) is the probability 
of assigning a galaxy that is at true redshift z to the 
photo-z bin between z p — a and b. We approximate 
the photo- z error to be Gaussian with bias 8z and rms 
&z = 020(1 + z ), an d the probability becomes 



T^p,i, 2p,<i *) = !Ki, z)/m z), 



B E 



I(a,b;z)- 



27r er. 



dz p exp 



(z. p — z — 8zY 



2(7? 



The normalization 1(0, 00; z) implies that galaxies with 
a negative photo-z have been excluded from n(z). Al- 
though uncertainties in the photo-z parameters 5z and 
a z have a l arge impact on the d ark energy cons t raints 
from WL (jHuterer et all 120061 : iMa et alJ l^OOl IZhanl 
2006), such systematic effects are not relevant to how 
WL probes dark energy. Hence, we fix 5z — and 
a z = 0.05(1 + z) per galaxy. In Section[5]we show results 
with photo- z and power spectrum systematics. 

Assuming that the observables, e.g., the shear map 
-f(0) in multipole space, and the likelihood of the pa- 
rameters of interest, q, near the fiducial model are both 
Gaussian, one can use the Fisher information matrix to 
estimate the errors of the parameters from the covari- 
ance, C, of the o bservables. In summary, the Fisher 
matrix is given by (|Teg mark 19971) 



v 2£+l Trr -!dC e x dC t 



(2) 



For BAO and WL tomography, we have 

where 6^ is the Kronecker delta function, 7 rms ~ 0.2 is 
the rms shear per galaxy for WL, and g rms = 1 for BAO. 
The Gaussian likelihood of the parameters is 

1 T 
-(q-q f ) F(q-q f ) 



C(q) cx If" 1 ! 1/2 exp 



where q^ corresponds to the fiducial model and F 1 
is the covariance of the parameters. The minimum 



marginalized la error of q a is a(q a ) — (F^ 1 )]!^ . In- 
dependent Fisher matrices are additive; a prior on q a , 
ap(q a ), can be introduced via F' aa — F aa + ap 2 (q a ). A 
Fisher matrix of the parameters q can be projected onto 
a new set of parameters p via 



F' 



(3) 



3. DISTANCE AND GROWTH FACTOR CONSTRAINTS 

The distance-redshift and growth factor-redshift rela- 
tions are the two most important diagnostics of the driv- 
ing mechanism of the accelerated cosmic expansion, be 
it dark energy or modified gravity. Various dark energy 
probes are essentially probes of these relations (and hence 
they are modified gravity probes as well) . Most models of 
dark energy and modified gravity do not modify the form 
of equation ([T]) . As such, reconstructing the distance and 
growth factor from the shear and galaxy power spectra is 
fairly model-indep endent, which is useful for distinguish- 
ing those models (ILue. Scoccimarro. fc Starkman] 120041: 
Sondl2005l: llshak. Upadhve. fc Sperge]||2006t iKnox et al l 
2006al Uain fc Zhand |2007t but cf. iLinderi 12001 
Bertschinger fc Zukinll2008l ). 

To reconstruct the distance and growth factor, one 
would set a series of distance and growth parameters cov- 
ering a range of redshifts, interpolate the continuous dis- 
tance and growth functions from a realization of these 
parameters, assign a probability of the observed data 
given the particular set of distance, growth, and other 
parameters, and repeat the process for enough realiza- 
tions to map the posterior distribution of all the param- 
eters. This is rather involved. However, if the purpose 
is to roughly estimate the errors without the actual ob- 
servational data, one often approximates the posterior 
distribution of the parameters as a multivariate Gaussian 
around its peak and applies the Fisher matrix analysis to 
infer the errors. We show here an example of estimating 
the errors of reconstructed distances and growth factors 
from BAO and WL, both separately and in combination. 

3.1. Parameters 

There is a minor ambiguity in the growth factor of 
the large-scale structure, and we choose the convention 
that G(z) = (1 + z)^ 1 in an Einstein-de Sitter universe 
but G(0) 7^ 1 in general. For distance parameters, it 
is convenient to use the comoving distance D(z) instead 
of the comoving angular diameter distance Da(z). The 
reason is that while D(z, z') = D(z') — D(z) holds all 
the time, Da(z, z') ^ Da(z') — Da(z) in the presence of 
curvature: 

K^^smlDiz^^K 1 / 2 ] K>{) 
D A {z,z') = { D(z,z') K = , 

\K\-^ 2 smh[D( Zl z')\K\^ 2 } K<0 

( 4 ) 

where K = —Qk{Ho/c) - Separately, the freedom of the 
Hubble parameter, H(z) = c[dD j dz)~ x , in the galaxy 
power spectra of equation (fTJ) could be unphysically re- 
stricted by the interpolation scheme. A conceptually 
simple remedy is to replace distance parameters Di with 
Hubble parameters Hi and then project Hi constraints 
onto Hq and Di constraints. 
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We assign 15 Hubble parameters Hi and 15 growth 
parameters Gi (i = . . . 14) at rcdshifts evenly spaced in 
log(l + z) from zq = to Zu = 5. The Hubble parameter 
H(z) and growth factor G(z) are then spline- interpolated 
from HiS and GiS. The comoving distance parameters 
Di are assigned at the same redshifts except that Dq = 
is replaced by Hq. Note that our distance parameters 
are in units of Mpc without h. For brevity, we only 
present the projected results without referring to Hi and 
the intermediate step. 

Given a cosmological model for the distance and 
growth factor, the parameter 9 S is a function of other 
explicit and implicit parameters that affect the distance 
and linear size of the sound horizon at the last scatter- 
ing surface. This means that, in the absence of the dis- 
tance and growth parameters, the numerical derivatives 
in equation @ with respect to each parameter must be 
taken with the implicit parameters h and fix perturbed 
in unison to keep other explicit parameters fixed. A con- 
sequence is that even though the physical matter density 
u> m , baryon density u>b, and 9 S do not directly affect the 
distance and growth factor, they do so indirectly through 
the implicit parameters. As one incorporates the dis- 
tance and growth parameters, the physical model that 
links the cosmological parameters to the distance and 
growth factor is no longer in effect, so information about 
the distance and growth factor must be carefully removed 
from the cosmological parameters. 

Parameters that affect the distance or growth factor 
without altering the normalization or the shape of the 
matter power spectrum may not be included in the fore- 
cast. Such parameters include the dark energy EOS pa- 
rameters and 9 S . The curvature parameter Qk has a 
negligible effect on the slope of the matter power spec- 
trum. It affects the angular diameter distance through 
its impact on the comoving distance, which is replaced 
by -DjS, and via equation which is a property of the 
Friedmann-Robertson- Walker (FRW) metric. As such, 
we keep 0^ for (and only for) calculating Da from D^s 
using equation Since the curvature parameter here 
has no role in the matter power spectrum, Di, or Gi, 
we label it as VL' k to distinguish from the real curvature 
parameter. 

Our approach is similar to that in iBernsteinl (|2006l ) 
where Da(z,z') is expressed in terms of Da(z), Da(z'), 
and, to first order, Q,' k . The constraint on Q! k will hold for 
any model that preserves the FRW metric and (the form 
of) equation JT]), whereas curvature constraints from ex- 
ploiting the full functional dependence of the angular 
diameter distance or luminosity distance on flk (e.g., 
ISpereel et al"1l200l are valid only for a particular cos- 
mological model. For this reason, measurem ents of f21 
are c onsidered pure metric tests for curvature (|Bernsteinl 
f2006h . 

The rest of the cosmological parameters do not affect 
the distance or growth factor. However, if the growth 
factor is already accounted for by the code that calcu- 
lates the matter power spectrum or trans fer function 
(e.g., CMBFAST. IZaldarriaga fe Seliakl 12000ft . one must 
use proper combinations of the code input parameters 
to obtain the derivatives for the Fisher matrix without 
altering the growth factor. 

In summary, our parameter set for forecasting distance 
and growth errors includes cu m , oJbi j Tj ? 3 



A% H , Di (i = 1...14), and G t (i = 0...14). For 
BAO, additional 15 linear galaxy bias parameters bi are 
included, which are assigned at the same redshifts as the 
growth parameters. In actual calculations, we use In Ho, 
In Di, lnGi, and lnfr^ as parameters, so their constraints 
are reported in fractional form. 

We assume a fiducial model for the linear galaxy bias 
b(z) = 1 + 0. 84,2, which is est i mated from the simula- 
tion results in Weinbe rg et al.l (|2004f ). The exact value 
of b{z) is not important for our purpose, though a higher 
bias does produce stronger signals (galaxy power spec- 
tra) and hence tighter parameter constraints. The de- 
pendence of the dark energy EOS e rror on the power 
spectrum amplitude can be found in[Zhan| (|2006ft . The 
nonlinear matter power spectrum is calculated using the 
iPeacock fe Dodda (|1996| ) fitting formula. A direct appli- 
cation of the fitting formula to the CDM power spectrum 
would cause a large shift of the BAO features. In addi- 
tion, it has difficulty processing power spectra that have 
an oscillating logarithmic slope. Thus, we first calcu- 
late a linear matter power spectrum with no BAO fea- 
tures that otherwise m atches the CDM power spectrum 
(|Eisenstein fe Hulll999f ). then take the ratio between the 
nonlinear and linear matter power spectra with no BAO 
features, and finally apply this ratio to the linear CDM 
power spectrum t o obtain the nonlinear CDM po wer 
spectrum (see also Eisenstcin et al. 2005; Zhan 2006). 

We modify cswab 3 (|Zhanl l2006t Izhan et al.N2008li to 
calculate the Fisher matrices for 10-bin (0 < z p < 3.5) 
WL and 20-bin (0.15 < z p < 3.5) BAO tomographic 
measurements over 20,000 square degrees. To reduce the 
impact of nonlinear evolution as well as ba ryonic effects 
on the small-scale matter power spectrum dWhitel 



Zhan fe Knoxll2004t lHagan et all 120051: Uing et alj 



2004; 



Rud d et al.1 12008). we limit the maximum multipole 



2006; 



to 2000 for WL and 3000 for BAO [the latter must also 
satisfy the condition Ag(£/DA', z) < 0.4]. At the low mul- 
tipole end, we exclude modes with £ < 40, which could be 
affected by dark ene rgy clustering (|Hu fc Scrantonll2004l ; 
ISong fc Knoxi r2004). Since Di, Gi, and 6j are not well 
constrained by the data at the low- and high-redshift 
ends, we impose extremely weak priors of ap(\nDi) = 
apQnGi) = apQnbi) = 100, i.e., 10,000% error, just to 
avoid vanishing diagonal elements in the Fisher matrices. 

3.2. Distance and Growth Factor Constraints 

How well can LSST measure the distance and growth 
factor given its statistical errors? Fig. [1] presents the es- 
timated errors on the distance (left panel) and growth 
(right panel) parameters from LSST WL (solid lines and 
open circles), BAO (dashed lines and open squares), and 
joint BAO and WL (dotted lines). The error of each 
distance (growth) parameter is marginalized over all the 
other parameters including the growth (distance) param- 
eters and other distance (growth) parameters. The lines 
show the results with weak independent Gaussian priors 
on cosmological parameters: o-p(lnuj m ) = o~p(hiuji,) — 
a P (SV k ) = a P (n s ) = a P (a s ) = 0.05, o~p (t) = 0.01, 
o-p(Yp) = 0.02, and cr P (lnAfj) = 0.1. These priors 
are rather weak compared to what can be achieved with 
Planck, except that fl' k is not constrained by CMB unless 
a particular cosmological model and appropriate priors 

3 Available at http://software.hzhan.net 
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Fig. 1. — Marginalized ltr errors on the comoving distance (left panel) and growth factor (right panel) parameters from LSST WL (solid 
line and open circles), BAO (dashed line and open squares), and joint WL and BAO (dotted line). The maximum multiple used for WL 
is 2000, and that for BAO is 3000 [with the additional requirement A$(£/ Da', z) < 0.4]. The growth parameters, Go • ■ ■ £ ) 14, are evenly 
spaced in log(l + z) between z = and 5, and the distance parameters, Di . . . D14, start at z\ = 0.14 (see text for details). The error of 
each distance (growth) parameter is marginalized over all the other parameters including growth (distance) parameters and other distance 
(growth) parameters. Lines represent results with cosmological parameters floating, i.e., marginalized with weak priors, and symbols are 
for results with the cosmological parameters fixed. 

as ^ ax . This alone could claim a factor of 3 increase in 
the error as ^ max decreases from 3000 to 1000. The small 
difference of the errors at low redshift is due to the fact 
that the highest multipole there is already limited by the 
condition Aj^/D^; z) < 0.4, a measure to exclude fully 
nonlinear modes. On the other hand, if we relax the con- 
ditions to Ag(£/DA; z) < 1 and £ max = 6000, the errors 
on the growth factor reduce to ~ 0.04 between z — 0.29 
and 3.1. 

It may seem quite optimistic to rely on the nonlinear 
power spectrum, even if restricted to mildly nonlinear 
regime, for determining the growth factor. However, this 
is not a serious concern, because the growth factor con- 
tributes little to BAO constraints on dark energy (see 
Section |4~T|) . Moreover, it has been demonstrated that 
there is no significant difference between the WL dark en- 
ergy constraints with a scale- free matter power spectrum 
or a nonlinear CDM power spectrum if the shot noise 
is negligible and that nonlinear evolution helps WL only 
because it boosts the signal on small scales where the sta- 
tistical uncertainty is limited by the shot noise (ZK06b). 
Nevertheless, Fig. Q] suggests that the galaxy power spec- 
tra could potentially constrain dark energy and alterna- 
tive gravity models with both distance and growth factor 
measurements as WL does if one understood the galaxy 
bias and nonlinea r evolution well through numerical 
simul a tions (e.g., ICen fc Ostrikerl l2000; Wei nberg et al.1 
120041: iHeitmann et all l2005h ~ semi-analytic model- 
ing (eXTTMafeFrvi [200fT iPeacock fc Smithl 120001 : 
Croton et al.1 l2007f) per turbative calculations (e.g., 
Jain fc Bertschingerlfl99l Ueong fc Komatsul f2006). and 



are assumed. In addition, they do not introduce extra 
correlations between the parameters. The errors on the 
galaxy bias parameters (not shown) trace the errors of 
BAO growth parameters in this case. The symbols as- 
sume that the cosmological parameters are known pre- 
cisely. A discussion on the sensitivity of these distance 
and growth factor constraints to individual cosmological 
parameters is given in Appendix IA"1 

It is not surprising that BAO places tighter con- 
straints on the distanc e parameters than WL (see also 
Seo fc Eisens"tehl l200l ISongl l2005t iKnox et all [2006al : 
Zhan fc Knoxl l2006al for results of either technique), 
but one may not expect that the former determines the 
growth parameters better than the latter at z > 0.7. It 
has been demonstrated that photo- z BAO measurements 
(without cross-bin power spectra) can constrain the com- 
bination of th e galaxy bias and grow th factor biGi to sev- 
eral percent (jZhan fc Knoxl [2006ah . but the degeneracy 
between hi and Gi could have left them undetermined 
individually. However, hi and Gi are not completely de- 
generate here; the linear galaxy bias is treated merely 
as a multiplicative factor in the galaxy power spectrum, 
whereas the growth factor (among others) determines 
both the amplitude and the shape of th e nonlinear mat- 
ter power spectrum (jPeacock fc DoddsllrOfll . Since the 
BAO technique is sensitive to the shape of the power 
spectrum, it is able to partially break the degeneracy 
and constrain both Gi and hi. This also means that the 
BAO constraints on the growth factor can be sensitive to 
the number of nonlinear modes included. For example, 
with £ max = 1000 and all the other parameters floating, 
the BAO error on the growth factor at z = 3.1 in the 
right panel of Fig. Q] increases from 0.064 to 0.21, though 
that at z — 1.2 increases by only 15%. The impact of the 
nonlinear information is compounded by the fact that in- 
cluding more modes, regardless whether they are linear 
or nonlinear, always improves the constraints on non- 
degenerate parameters. If the trace term in equation j2]) 
is roughly scale-independent, the Fisher matrix will scale 



direct measurements fr om higher - order statistics of 
the galaxy distrib utio n fFrvl Il994t iMo. Jing. fc White! 
19971 IVerde et al.ll2002l; iGaztanaga fc Scoccimarroll2005l ; 
Sefusatti fc Komatsufl2007l ). 

The joint BAO and WL analysis reduces the errors on 
the distance and growth factor parameters. Although 
the improvement on the marginalized distance errors is 
moderate, the resulting constraints on the dark energy 
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EOS parameters are much tighter than those from either 
technique alone (see Section |4~2|) . The marked improve- 
ment on the growth factor parameters is due to the fact 
that the shear power spectra, galaxy-shear power spec- 
tra, and galaxy power spectra have different dependen- 
cies on the galaxy bias (6°, b, and b 2 , respectively). This 
further lifts the degeneracy between Gi and bi, allowing 
significantly better determinations of Gi . 

3.3. Metric Test of Curvature 

The mean curvature of the universe is of great theoret- 
ical interest. It is estimated that future large-scale BAO, 
SN, and WL surveys can constrain ilk (or fifcft 2 ) to bet- 
ter than 10~ 3 with the assumption of matter dominance 
at z > 2 and precise independent distance me asurements 
at z > 2 and at recombination (Knox 2006) or with a 
specific dark energ y EOS w(z) = wq + w a z(l + z)~ x 
(|Knox et all l2006ri IZhanl [2006). However, if one as- 
sumes only the Robertson- Walker metric without invok- 
ing the dependence of the comoving distance on cos- 
mology, then the pure metric constraint on curvature 
from a simple combination of BAO and WL becomes 
tr(n' fc ) ~ 0.04/ s ^ r /2 (a z0 /0.04) 1 / 2 (lBernsteinl l2006l). 

Our result for Q' k from LSST WL or BAO alone is 
not meaningful, but the joint analysis of the two leads 
to a(£l' k ) = 0.015 with the aforementioned weak priors 
except that no prior is applied to Q' k . The error on Q' k 
changes only mildly with the number of parameters used, 
e.g., a(Q' k ) = 0.017 for 100 ds and 99 As, even though 
errors of (the eigenmodes of) the distance and growth 
factor scale roughly as the square root of the number 
of parameters. This is expected, because the mapping 
between the underlying comoving distance and the ob- 
servable angular diameter distance through Q' k does not 
depend on the number of parameters. 

Planck will provide an accurate measurement of the 
angular diameter distance to the last scattering surface, 
and it can reduce the error moderately from a(fi' k ) = 
0.015 to 0.013. Separately, the constraint weakens to 
= 0.017 with conservative estimates of system- 
atics in photo-zs and power spectra expected for LSST 
(see Section [5] for more details). This is better than the 
forecast derived from the shear po wer spectra and galaxy 
power spectra in lBernsteinl ()2006l ). because we include in 
our analysis more information: the galaxy-shear power 
spectra. 

4. PROBING DARK ENERGY WITH DISTANCE AND 
GROWTH FACTOR 

With the distance and growth factor constraints from 
the last section, one can now answer the questions about 
the roles of the distance and growth factor in probing 
dark energy and the difference between the distance mea- 
surements from BAO and WL. 

To assess the utility of distance and growth fac- 
tor in constraining dark energy, we adopt a generic 
parametrization of the dark energy EOS w(z) — wq + 
w a z{l + z)^ 1 and project the distance and growth factor 
errors onto wq and w a . We use the product of la errors, 
a(w p ) x <j(wa), to gauge dark energy constraints, where 
the er ror a(w r> ) is equal to that on wp with w a b eing held 
fixed (|Huterer fc Turnerll2001l : lHu & Jamll2004h . This er- 
ror product (EP) is proportional to the area of the error 
ellipse in the wo—w a plane. The figure of merit in the 
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Fig. 2. — WL la error contours of the dark energy EOS param- 
eters Wd and w g , which artificially split the EOS into a component 
that affects the distance only (w^) and another that affects the 
growth factor only (w g ). If all the other cosmological parameters 
are allowed to float with weak priors (dot-dashed line), the con- 
straint on w,i from Di will be much stronger than that on w g from 
Gi, whereas if all the other cosmological parameters are fixed (solid 
line), the constraints on w^ and w g will become comparable. 

Dark Energy Task Force report (Al brecht et al.ll2006l) is 
the reciprocal of the product of 2a errors on w p and w a - 

4.1. Roles of Distance and Growth Factor 

The relative power of the distance and growth factor 
in constraining the dark energy EOS has been studied in 
a number of ways. For example, SB05 and ZHS05 split 
the dark energy EOS into two constant components: Wd 
that affects only the distance and w g that affects only 
the growth factor. They find that the constraints a(wd) 
from the distance and a(w a ) from the grow th factor are 
comparable. However. iKnox et al.1 (j2006aT ). by project- 
ing estimated distance and growth factor errors onto a 
constant EOS w, mention that the constraint on w is 
due mostly to the distance. ZK06b explores this subject 
in yet another way; they compare the EPs with cosmo- 
logical information removed from either the distance or 
the growth factor. Their approach is equivalent to fix- 
ing Di (d) and then projecting the errors of Gi (Di) 
onto cosmological parameter space, though they bypass 
the intermediate stage of estimating Di and Gi errors. 
ZK06b shows that WL distances are much more power- 
ful than WL growth factors in constraining wq and w a 
if other cosmological parameters are marginalized with 
weak priors and that they are comparable if the other 
cosmological parameters are fixed. 

Because of the different methods employed, it is not 
straightforward to synthesize the above mentioned re- 
su lts. Our mu l ti-stage d approach is compatible with that 
in lKnox et al.l (|2006af ). and we obtain fully consistent re- 
sults with those in ZK06b when their method is applied 
to the Di and Gi errors in Section l3~2l Hence, we only 
need to explore the results for the method in SB05 and 
ZHS05. 

We project the WL Fisher matrix in Section l3~2l onto 
(w d , w g , u m , w b , 9 Sl fifc, r, Y p , n s , a s , A 2 R ) using 
equation (J3]) with dDi/dw g = dGi/dwd = 0. The re- 
sulting la error contours of Wd and w g are shown in 
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TABLE 1 

Marginalized la Errors on Dark Energy 
Parameters with All the Other 
cosmological parameters floating 



Probe 


Di 


G t 


Wo 


w a 


w p 


EP" 1 


WL 


Y 


Y 


0.078 


0.22 


0.023 


200 




Y 


M 


0.23 


0.71 


0.054 


26 




M 


Y 


1.4 


2.8 


0.23 


1.6 


BAO 


Y 


Y 


0.33 


0.81 


0.044 


28 




Y 


M 


0.34 


0.85 


0.045 


26 




M 


Y 


3.1 


5.5 


0.87 


0.21 


Note. 




"Y" 


indicates that the set of 


par am- 



eters are projected onto cosmological parameter 
space, and "M" indicates that they are marginal- 
ized. Weak priors are applied to cosmological pa- 
rameters except wq and w a , as discussed in Section 
[3721 For both BAO and WL, the distance is far 
more powerful than the growth factor in constrain- 
ing wo and w a - 

Fig. |U with a(w d ) = 0.036 and <r(w g ) = 0.15 when all 
the other cosmological parameters are marginalized with 
weak priors (dot-dashed line) and a(wd) = 0.0080 and 
<j{w g ) = 0.0071 when all the other cosmological parame- 
ters are fixed (solid line) . The latter result is roughly con- 
sistent with <j(wd)/<j(wg) ~ 1/2 from a single-bin shear 
power spectrum in SB05 and a(utd)/cr(w g ) ~ 1 in ZHS05, 
though the orientation of the error contour in ZHS05 is 
nearly orthogonal to that in Fig. [2] and that in SB05. To 
examine if the differences of the assumed surveys could 
explain this discrepancy, we re-calculate the constraints 
with the survey parameters used in SB05 (f s ky = 0.024, 
a — (3 — 2, z* = 1.13, n to t = lOOarcmin -2 , l max = 10 4 , 
no tomography, & no redshift error). The resulting er- 
rors, <j(wd) — 0.027 and cr(w g ) = 0.032 with all the other 
parameters fixed, and the orientation of the error contour 
are in good agreement with those in SB05 except that our 
<j(w g ) is ~ 40% smaller. We are not able to reproduce 
Fig. 3 of ZHS05 (/ sky = 0.1, a = 2, f3 = 1.5, z* = 1, 
"■tot = lOOarcmin -2 , £ max = 200, 20-bin tomography, 
& <j z q = 0.05); our errors on Wd and w g remain pos- 
itively correlated with comparable size [<j{wd) = 0.051 
& cr{wg) — 0.081] when all the other parameters are 
fixed and anti-correlated with much smaller error on uid 
[a(wd) — 0.25 & cr(w g ) = 1.2] when all the other param- 
eters are marginalized with weak priors. We note that 
even though Planck will place very tight constraints on 
CMB-related parameters, it does not have the same effect 
as fixing those parameters, e.g., we obtain a(wd) — 0.022 
and a(w g ) = 0.099 for LSST with Planck priors. There- 
fore, the distance should be at least several times more 
powerful than the growth factor in constraining the dark 
energy EOS for WL tomography. 

The parameter splitting technique has been applied to 
CMB, galaxy clustering, SN, and W L data to check th e 
consistency of dark energy models (|Wang et al.l l2007h . 
The results do show that the error on w g is several times 
larger than that on Wd, even without the galaxy clus- 
tering and SN data that provide good distance measure- 
ments at low redshift. 

Another approach to isolate the roles played by the 
distance and growth factor is to separately marginalize 
distance or growth parameters and evaluate the dark en- 
ergy constraints from the rest. Tables [T] (with all the 



TABLE 2 

Marginalized 1<t Errors on Dark Energy 
Parameters with All the Other 
Cosmological Parameters Fixed 



Probe 


Di 


Gi 


w 


w a 


w p 


EP" 1 


WL 


Y 


Y 


0.062 


0.15 


0.0069 


990 




Y 


M 


0.12 


0.32 


0.029 


110 




M 


Y 


0.84 


1.8 


0.036 


15 


BAO 


Y 


Y 


0.19 


0.39 


0.022 


120 




Y 


M 


0.20 


0.42 


0.022 


110 




M 


Y 


1.3 


2.3 


0.072 


6.2 


Note. 




Same 


as Tabic [TJ but 


with all the other 



cosmological parameters fixed. 



other cosmological parameters floating) and [2] (with all 
the other cosmological parameters fixed) present such 
tests. There is no doubt that the growth factor is im- 
portant for probing dark energy with WL tomography, 
as, in both tables, the EP increases by a factor of ~ 8 if 
the growth parameters are marginalized. However, given 
the increase of EP by a factor of ~ 100 when the distance 
parameters are marginalized, the growth factor is clearly 
much less constraining on the dark energy EOS than the 
distance. 

Tables [1] and [2] include BAO results as well. We see in 
both tables that the dark energy constraints are weak- 
ened by less than 10% when the growth parameters are 
marginalized but more than an order of magnitude when 
the distance parameters are marginalized. This demon- 
strates that the growth factor plays almost no role in 
constraining wq and w a with photo-z BAO, even if it 
can be measured from the shape of the nonlinear power 
spectrum (Fig.[T]). 

Comparing the separate BAO and WL constraints 
from Di in Tables Q] and [2j we find that for our fidu- 
cial survey, BAO and WL achieve similar EPs (though 
different errors on wo and w a individually) in the absence 
of the growth information (and systematics) . More inter- 
estingly, WL Gi provides more information on the dark 
energy EOS than BAO Gi, despite that the errors of the 
latter are smaller than those of the former at z > 0.7. 
The reason is that the growth factor is less sensitive to 
wq and w a at higher redshift and that the low-redshift 
growth factors, which WL measures better than BAO, 
are helpful in breaking the degener acy betwee n and 
w a (as are high- redshift distances) (|Zhari 2006; ZK06b). 

To gain a further understanding of the roles of the 
distance and growth factor in constraining dark energy, 
we show in Fig. [3] marginalized la contours of AlnG;- 
A In Gj (lower triangle) and A In Di~A In Dj (upper tri- 
angle) for WL (solid lines) and BAO (shaded areas). 
The line segments in the figure indicate the amount of 
change, AlnG^ and AlnDi, due to small perturbations 
Aw = ±0.1 (dashed lines) and Ail k = 0.005 (dotted 
lines). The sensitivities of In Gi and In Di to w a (not 
shown) are roughly 1/2 and 1/3 of those to w , respec- 
tively. From these line segments, one can see that the 
sensitivity of the growth factor to fifc and that of the 
distance to wq (and w a as well) decrease quickly with 
increasing redshift. The sensitivity of the growth factor 
to wo and w a has a shallow peak around z ~ 0.6. 

Fig. [3] illustrates a major advantage of WL growth 
factors over BAO growth factors, if ever measured, in 
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Fig. 3. — Comparison vs. redshift of WL (solid lines) and BAO (shaded areas) constraints on the distance (upper triangle) and growth 
(lower triangle) parameters and sensitivities of InD; and InGi to u>o and Q^. The error ellipses are lcr and marginalized over all other 
parameters including irrelevant DiS and GiS. The line segments correspond to the changes, A In Gi— AlnGj and A In Di— A In D 3 ;, due to 
small perturbations Amio = ±0.1 (dashed lines) and Afife = 0.005 (dotted lines). Note that AlnG^-AlnGj are enlarged 10 times for Awrj 
and 5 times for Af2fe. Though not shown, the amount of change, AlnGi (AlnD^), with Aw a = ±0.1 is roughly half (one third) of that 
with Auiq = ±0.1. The WL growth constraints are tighter at lower redshift where the parameter sensitivities are higher. Moreover, the 
WL growth error ellipses are often oriented in the most sensitive direction to the parameters, so that the WL growth constraints provide 
more dark energy information than the BAO ones (Tables [T] & [2j , even when the marginalized BAO growth errors are smaller (Fig. [TJ- 



constraining the dark energy EOS: WL determines the 
growth factor more accurately than BAO at low redshifts 
where parameter sensitivities are higher. Another promi- 
nent difference between the BAO and WL constraints on 
the growth factor is that the WL A In Gi— A In Gj con- 
tours are often oriented in the most sensitive direction 
to the parameters, i.e., orthogonal to AlnG^-AlnGj 
caused by Awq or AQ&, whereas the BAO ones are gen- 
erally in the least sensitive direction. Therefore, one can 
still extract more information from the WL Gi even when 
the marginalized errors of BAO Gi are smaller. Such an 
advantage is seen up to z ~ 2 beyond which the WL 
contours start to enclose the BAO contours completely, 



but the parameter sensitivities also become fairly low at 
z > 2. 

4.2. Insights from Distance Eigenmodes 

Fig. Q] and Tables [1] and [5] pose another puzzle: while 
the BAO distances are more accurate than the WL dis- 
tances, the errors on wq and w a from the BAO distances 
are actually larger than those from the WL ones. ZK06b 
postulates that the correlation of the WL distance errors 
might hold the answer, but, unlike the growth factor er- 
ror contours, the WL and BAO distance error contours 
in Fig. [3] are generally in the same direction with the 
former enclosing the latter. To resolve this puzzle, one 
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has to realize that the contours in Fig. [3] are marginal- 
ized over all the other parameters including irrelevant 
distance and growth parameters, and useful information 
may have been lost in the process. To examine such in- 
formation, an cigenmode analysis is needed. 

Let Cd be a distance covariance (sub-)matrix from 
Section 13.21 By definition, it is marginalized over all 
the other parameters with their priors. The distance 
eigenmodes Ej = (E n ,E I2 , . . . , E ItNp _ l ) T satisfy 

Cy>Ej = Xj E j I 



1. 



1, 



(5) 

where A/ is the eigenvalue corresponding to Ej, N p — 1 
is the number of distance parameters (remember that 
Dq = is replaced by Hq), and capitalized indices are 
used to avoid confusion with the indices of the distance 
parameters. For convenience, we sort Ei by their eigen- 
values in ascending order. Since the eigenmodes are or- 
thonormal to each other, their covariance matrix, A, is 
diagonal with the variance <r 2 (Ei) = A/, i.e., 

Ajj = SfjXj. (6) 

One can interpret these Ej as the modes of departure of 
the distances from the fiducial model and A/ as a measure 
of how well such deviations can be determined, regardless 
the driving mechanism. 



With equations ((5]) and (' 
ized distance Fisher matrix 



y), we project the marginal- 
into WQ-Wa space 



F a (3 — 



Np-1 

E 

- E 



d In Di 

dp a 

d In Di 

dp a 



din Dj 

dp/3 



EnXr 1 E 



dbxDj 



(7) 



with p — (wo,w a ) T . We see that the contributions of 
the distance eigenmodes to the final Fisher matrix are 
separable 



Fa/3 — 



£la£ip 

7=1 



A, 



(8) 



where Ej a = Y). En d\n.Dj/dp a . In other words, the 
WQ-w a Fisher matrix consists of statistically independent 
contributions from the distance eigenmodes. 

The number of distance parameters, N p — 1, does not 
have much impact on the projected dark energy con- 
straints, but it can assert artificial restrictions on these 
modes. Thus, before drawing conclusions from the dis- 
tance eigenmodes, one must ensure that the interpolation 
has sufficient degrees of freedom to reflect the intrinsic 
properties of the distance eigenmodes of the dark energy 
probes, though having too many distance parameters will 
leave the results prone to numerical errors. Since the in- 
trinsic eigenmodes do not vary with N p , one can increase 
N p until the shapes of the eigenmodes converge. In addi- 
tion, the variance a 2 (Ei) should grow in proportion with 
the degrees of freedom, N p — 1, once the convergence is 
reached. By the same token, we change the priors on 

1/2 

hiDi, InGi, and \nbi to 100N P , though these priors 
are already too weak to make any difference. 

Fig. H] shows the errors of the distance eigenmodes 
scaled by (N p — l)~ x / 2 for WL (solid line and filled sym- 
bols) and BAO (dotted line and open symbols) with 
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Fig. 4. — Comparison of the errors of the distance eigenmodes, 
cr(Ej), from WL (solid line and filled symbols) and from BAO (dot- 
ted line and open symbols) . All the other parameters are marginal- 
ized, and the errors are scaled by the square root of the number of 
distance parameters, Np — 1, with Np = 15 (circles), 50 (squares), 
and 100 (lines). The WL distances have several eigenmodes that 
are better determined than all the BAO distance eigenmodes. The 



alignment between the N p 



50 and Np 



100 results suggests 



that the best determined modes have converged at N p > 50. 

N p = 15 (circles), 50 (squares), and 100 (lines). The 
alignment between the N p = 50 and N p = 100 eigenmode 
errors suggests that N p > 50 is sufficient for the best de- 
termined modes to converge, which is confirmed upon 
inspection of the modes. Thereafter, we set N p = 100. 

An interesting result in Fig.[4]is that some WL distance 
eigenmodes are better determined than all the BAO dis- 
tance eigenmodes, even though the marginalized errors 
of individual WL distance parameters are all larger than 
those of the BAO distance parameters in Fig. [T] This 
is predicted in ZK06b to be an advantage of the WL 
distance measurements, albeit that the best determined 
modes do not necessarily provide the most information 
on dark energy. 

Fig. [5] shows a small piece of the distance covariance 
matrix of WL (left panel) and that of BAO (right panel) . 
The covariance matrices are normalized 



(C D )i 



(C D ) i 



so that |(Cb),j| < 1. Because the lensing kernel is much 
broader than the galaxy power spectrum kernel, we ex- 
pect the WL distance covariance to vary more slowly 
than the BAO distance covariance, which is indeed seen 
in Fig. [5J The features of the distance covariance matri- 
ces in Fig. [5] are typical over the whole redshift range, 
with some broadening at both low-z and high-2: ends, 
though the values of the individual covariance elements 
do vary with the details of the calculations such as the 
interpolation scheme and tomographic binning. 

The characteristic difference between the distance co- 
variance matrices of the two techniques must hold the an- 
swer to the paradox brought up at the beginning of this 
subsection: although BAO determines the distance more 
accurately than WL (the WL distance error contours also 
enclose the BAO ones), the EPs of each technique, when 
the growth parameters are marginalized, are the same 
with WL having smaller marginalized errors on wq and 
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BAO (C D )ij 




Fig. 5. — Structure of the distance covariance matrix from WL (left panel) and that from BAO (right panel). The difference between the 
two covariances reflects the complementarity of BAO and WL distances. Elements of the covariance matrices are scaled by the geometric 
mean of corresponding diagonal elements. Only 16 distance parameters between z = 1.14 and 1.81 are presented here. The covariance 
matrices have been smoothly interpolated to show the underlying structure, and the discrete matrix elements are identifiable with the peaks 
and troughs if present. Note that the values of the matrix elements change rather slowly along the diagonal direction, so the apparent 
oscillations along that direction are not real. They are the result of isotropic rendering by the graphics software for a grid of data that vary 
more rapidly in the direction perpendicular to the diagonal than in the direction of the diagonal. 
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Fig. 6. — The best determined distance eigenmodes (left panels) and the ones that best constrain wo and w a (right panels) from WL 
(upper panels) and BAO (lower panels). The modes are sorted by their error in an ascending order. All the other parameters are floating, 
except that one of the best determined BAO distance eigenmodes (dashed line) is obtained with all the cosmological parameters fixed. The 
best determined WL distance eigenmodes do not change much whether the cosmological parameters are fixed or not. While WL is best at 
measuring distance modes that are insensitive to the absolute distance normalization (upper left panel), its constraints on wo and w a are 
derived largely from modes that are sensitive to the distance normalization (upper right panel). 



w a . Moreover, this difference also suggests that BAO 
and WL distance measurements are complementary. 

We present several eigenmodes of BAO and WL dis- 
tances in Fig. [5] While not shown, the best deter- 
mined WL modes (upper left panel) do not change much 
whether the cosmological parameters are fixed or not. 
These modes have close to zero mean, i.e., | J2i E n\ ^ 1> 
and, hence, are not sensitive to the absolute normaliza- 
tion of the distance. To see this, assume that all the dis- 
tances are subject to the same multiplicative factor 1 + / 
(| /| <C 1), so that the fractional errors A In Di = f. After 
projecting the distance errors on the mode E /, one gets 
the amplitude £\ AlnA-E/i = fY,i E Hi Le -> \T,i E n\ 



gauges the sensitivity of an eigenmode to a change in 
the distance normalization. Thus the best determined 
WL modes show that WL is best at measuring the rel- 
ative distance, which agrees with the expectation that 
WL tomography measures distance ratios from the lens- 
ing kernel. 

The best determined BAO distance eigenmodes (lower 
left panel in Fig. [6]) display more sensitivity to the ab- 
solute distance than the best determined WL distance 
eigenmodes, especially when the standard ruler of BAOs 
is precisely known by fixing the cosmological parameters. 
The best determined BAO modes peak at higher redshift 
than the WL ones because of the following reasons. First, 
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Fig. 7. — Decomposition of the wo and w a constraints from the smallest-EP distance eigenmodes (broken lines) for WL (left panel) and 
BAO (right panel). The WL distance eigenmodes are fairly complementary to each other, so that only 5 modes are needed to account for 
nearly all the dark energy constraints of the 99 distance eigenmodes (solid line). The BAO modes generally produce narrower wo— w a error 
contours than WL ones, but their orientations are similar. This makes the BAO modes less effective in reducing the marginalized errors on 
wo and w a individually, even though the error contours from all the distance modes (solid lines) of LSST BAO and WL have similar areas 
(or EPs, see Table[l]l. BAO and WL are highly complementary to each other, and the joint analysis (shaded area in both panels) reduces 
the EP by a factor of ~ 9 with respect to the EP of either technique alone. We have applied the priors ap(wo) = 1 and op(w a ) = 2 to 
individual modes and marginalized over all the other parameters with weak priors. The indices of the WL smallest-EP modes are I = 
7, 5, 17, 12, 8, 4, 14, 15, 23, in an order of increasing error ellipse area (or EP), and those of BAO modes are / = 29, 26, 27, 39, 31, 35, 38, 32, 23. 
The fact that the best determined distance eigenmodes from both BAO and WL are not the most sensitive modes to uiq and w a means that 
the two dark energy probes can potentially explore a larger dark energy parameter space and afford some redundancy to control certain 
systematic uncertainties. 



the severe truncations of multipoles of the galaxy power 
spectra at low redshift, e.g., f max = 732 at z p = 1.05 
vs. £ max = 3000 at z p — 2.09, increase the low-z BAO 
distance errors. Second, given the same source galaxy 
distribution, the lensing kernel peaks half-way between 
the observer and the source, whereas the galaxy power 
spectrum kernel peaks at the source. Finally, at higher 
redshift, the BAO features move to higher multipoles, 
where the smaller cosmic variance leads to better mea- 
surement of the distance. 

The best determined distance eigenmodes are not nec- 
essarily the ones that contribute the most to the con- 
straints on the dark energy EOS, because they may not 
be sensitive to a particular deviation of the distance- 
rcdsh ift relation caused by dark energy. Phrased another 
wav. lde Putter fc Linderl (|2007f) point out that the eigen- 
mode analysis gives the expected noise of the measure- 
ments, i.e., the error of the eigenmode, but the signal, 
i.e., the amplitude of the mode, and the signal-to- noise 
ratio are no t known unless com parison models are as- 
sumed (e.g.. iBarnard et al"1 l2008) or until measurements 
are made. 

For the wo-w a parametrization of the dark energy 
EOS, the comparison models populate the WQ-w a plane 
with the fiducial model (cosmological constant) at wo = 
— 1 and w a = 0. The contribution of each distance eigen- 
mode to the WQ—w a Fisher matrix is statistically indepen- 
dent of each other, so that we can evaluate the usefulness 
of an eigenmode by its constraint on wq and w a . Since a 
single mode cannot constrain both wq and w a , we apply 
the priors <jp(wq) — 1 and op(w a ) = 2 in the projection 
of a(Ei) onto wo and w a via equation ([8]). 

For a broad class of dark energy models, their effect 
on the distance is concentrated at low redshift without 
rapid oscillations. This means that a change of sign in 



En at low redshift will result in a reduction of the sen- 
sitivity to the dark energy EOS. Therefore, dark-energy- 
cfHcicnt modes should be roughly those with a large value 
of | J^i En\ at low redshift. The right column of Fig. [6] 
shows the distance eigenmodes with the smallest EPs 
for WL (upper panel) and BAO (lower panel). As ex- 
pected, they all have a major peak at z < 1. The BAO 
modes (lower right panel) exhibit oscillatory behavior at 
z > 0.6 because they are orthogonalized to other better 
determined modes that vary smoothly. 

The fact that the best determined distance eigen- 
modes (with other parameters floating) from both WL 
and BAO are not the most sensitive modes to wq and 
w a means that the two dark energy probes can poten- 
tially explore a larger dark energy parameter spa ce (e.g., 
lAlbrecht fc Bernstein! [20071 : IBarnard et ai1l2008f) and af- 
ford some redundancy to control certain systematic un- 
certainties. We also note that although cr{Ej) and EP 
do not have one-to-one correspondence, contributions to 
the dark energy constraints naturally come mostly from 
the best determined modes. For WL, the dark-energy- 
sensitive modes are confined within / < 20, whereas for 
BAO, they are spread out to / < 40. This is reflected in 
the indices of the smallest-EP modes in Fig. [6] as well. 

There has been a notion that WL tomography con- 
strains dark energy with relative distances or distance 
ratios from the lensing kernel. It is true that the best 
determined WL distance eigenmodes are not sensitive to 
the normalization of the distance, but the modes that 
contribute the most to the constraints on wq and w a 
are more sensitive to the normalization of the distance. 
For example, the average value of | En\ of the 5 best 
determined WL distance eigenmodes is 0.0094, whereas 
that of the 5 smallest-EP modes is 1.7 (for BAO, the val- 
ues are 0.033 and 1.4, respectively). Since the above re- 
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suits are marginalized over the growth parameters, which 
weakens WL's sensitivity to the absolute distance, it is 
reasonable to attribute the WL constrains on dark energy 
to its measurements of absolute distances. 

To see why the WL distances give smaller errors on wq 
and w a , even though they are less accurately determined 
than the BAO distances, we show in Fig. [7] marginalized 
la WQ-Wa error contours of the eigenmodes that have the 
smallest EPs of each technique. The priors ap(wo) = 1 
and <Jp(w a ) = 2 are applied, and all the other parameters 
are marginalized with weak priors. The solid error con- 
tour in each panel shows the constraints of all 99 distance 
eigenmodes for each technique. Since the dark energy 
constraint from each eigenmode is independent of those 
from other eigenmodes, one can see immediately that 5 
WL distance eigenmodes can account for nearly all the 
wo and w a constraints of WL distances. This echoes the 
conclusion from Fig. [6] that WL tomography has the po- 
tential of constraining more than wq and w a with all the 
well determined modes. Furthermore, the error ellipses 
of these 5 eigenmodes are oriented differently, so that 
they are effective in reducing the marginalized errors on 
wq and w a . The BAO error ellipses are narrower than 
the WL ones in general, but they are aligned in nearly 
the same direction. As such, combinations of the photo- 2 
BAO distance eigenmodes are not effective in reducing 
the marginalized errors of wq and w a , and many more 
modes are needed to match the combined EP with that 
of WL (Tables [Hand [I]). Because BAO and WL have dif- 
ferent parameter degeneracy directions and because the 
galaxy-shear cross power spectra provide additional in- 
formation, the two techniques are highly complementary 
to each other. Indeed, the joint BAO and WL distance 
constraints on wq and w a (shaded area) is roughly 9 times 
tighter tha n those from either technique alone (see also 
[Zhan|[200l) . 

5. CONCLUSIONS 

We have disentangled the roles of the distance and 
growth factor in constraining the dark energy EOS and 
given a detailed comparison between the distance mea- 
surements from BAO and WL as well as their constraints 
on wq and w a . We find that if the growth parameters 
are marginalized without being projected onto wq and 
w a constraints, the EP of LSST WL will increase by a 
factor of ~ 8. However, if the distance parameters are 
marginalized instead, the degradation will be two orders 
of magnitude. This clearly shows that the growth factor 
is important to WL, but the distance is far more impor- 
tant. In contrast, the BAO EP increases by only 10%, 
when the growth parameters are marginalized. It is also 
interesting that LSST BAO and WL achieve similar EPs 
when the growth parameters are marginalized. 

The reconstruction of a continuous function depends 
on how the function is represented by discrete param- 
eters. To explore the intrinsic properties of the recon- 
structed distances from BAO and WL, we assign a large 
number of distance and growth parameters and examine 
various aspects of the distance eigenmodes. This exercise 
leads to the finding that the WL distance eigenmodes are 
more complementary to each other than the BAO ones. 
As such, 5 WL eigenmodes can provide most of the LSST 
WL distance constraints on wo and w a , whereas many 
more are needed for BAO. A more useful result is that 
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Fig. 8. — Same as Fig.[T] but for distance (solid line and open tri- 
angles) and growth factor (dotted line and open circles) constraints 
from the joint analysis of LSST BAO and WL with a conserva- 
tive level of systematic uncertainties in the photo-z error distribu- 
tion and additive and multiplicative errors in the shear and galaxy 
power spectra (see text for more details). The lines represent the 
results with other parameters floating, and the symbols are with 
CMB priors from Planck. The joint constraints on distance are 
relatively insensitive to the assumed systematics. 

both LSST BAO and WL have some well determined dis- 
tance eigenmodes that are not very sensitive to wq and 
w a . These modes may be able to constrain additional 
dark energy parameters and some systematic uncertain- 
ties. 

The insight gained from the results in Sections [3] and 
[4] in the absence of various systematics is valuable for 
understanding the techniques. It is also of great interest 
for future surveys to have more realistic error estimates 
that include those systematics. To this end, we present in 
Fig. [D constraints on the distance (solid line and open tri- 
angles) and growth factor (dotted line and open circles) 
from the joint analysis of LSST BAO and WL with a con- 
servative level of systematic uncertainties in the photo-z 
error distribution and additive and multiplicative errors 
in the shear and galaxy power spectra. The lines repre- 
sent the results with other parameters floating, and the 
symbols are with CMB priors from Planck. By com- 
paring with Fig. [TJ one sees that the joint constraints 
on distance are relatively insensitive to the assumed sys- 
tematics whereas those on the growth factor are affected 
more. 

The results in Fig. [8] are marginalized over additional 
30 photo-z bias parameters 5zi, 30 photo-z rms param- 
eters <J z ,i, 10 shear additive error parameters AJ, 10 
shear multiplicative error parameters fj, and 20 galaxy 
additive noise parameters Af. We extend the addi- 
tive and mult i plicat ive shear power spectrum errors in 
iHuterer et al.l (|2006( ) to include the galaxy power spec- 
trum errors: 



(cf y )« = (i + W + fyfY)P$ r M + 



>XY 



cK y2 - 
"ij rms 



X A X \Y 



J_ 



where p x determines how strongly the additive errors 
of two different bins are correlated, and T] X and £ x ac- 
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count for the scale dependence of the additive errors. 
Note that the multiplicative error of galaxy number den- 
sity is degenerate with the galaxy clustering bias and 
is hence absorbed by hi. Below the levels of system- 
atics future surveys aim to achieve, the most impor- 
tant aspect of the (shear) additive error is its ampli- 
tude (|Huterer et all |2006[ ), so we simply fix p x = 1 
and r\ x = 0. For more comprehe nsive accounts of th e 
above syst e matic uncertainties, seefHuterer et al.l (|2006l ); 
I.Tain et all (12001 : iMa et al.l piM ): iZhanl (12006T ). 

We have applied priors ap(Sz) = 2~ 1 / 2 op{(j z ) = 
0.01(1 + 2:) per photo-z parameter interval (Az = 0.17) in 
Fig. El For Gaussian photo-z errors with a z — 0.05(1 + 
z), these priors correspond to a calibration requirement 
of 25 spectra per redshift interval of 0.17. The constraints 
do not change appreciably even if one tightens the pri- 
ors to 0.003(1 + z), because the photo-z sample-photo-z 
sample cross correlations ca n calibrate the uncertainties 
in Sz and a z to 10 -3 level (|Zhartl2006T ). Similar results 
are obtained per redshift bin with s pectroscopic sa mple- 
photo-z sample cross correla tions (iNewmanl 12 008). For 
the shear multiplicative error. lMassev et al.l (|2007l ) shows 
that current methods consistently achieve better than 2% 
precision. Allowing for another 5 to 10 years of develop- 
ment, we are hopeful that the multiplicative errors can 
be controlled to half a p ercent, i.e., ap(.f7) = 0.005. For 
the shear additive error, iPaulin-H cnrikss on et al.l (|2008l ) 
finds for a 15,000 deg 2 ground survey that it will be 
well under control over relevant scales if the point spread 
function (PSF) is calibrated over 100 stars. LSST will be 
able to use on average 2-3 stars per square arcminute 
to calibrate the PSF in each exposure, so the shear ad- 
ditive error will not have a significant impact on scales 
larger than ~ 5 arcminute (or I < 2000) even for a sin- 
gle ex posure. This is supported by a study with Subaru 
data (|Wittmanl 120051 ) . which shows that the correlation 
of the residual stellar shear after the PSF correction using 
only 0.9 star per square arcminute (out of ~ 8 arcmin -2 
available) is well below the cosmic shear signal at ar- 
cminute scale for a single 10-second exposure and that 
it is roughly inversely proportional to the number of ex- 



posures. In addition to deep exposures in m, g, and y, 
LSST plans to take 400 exposures per sky field in each 
of its riz filters (these will be used for WL as well as 
photo-z). Thus it is reasonable to project that the ad- 
ditive shear error will be sub-dominant to the statistical 
errors on relevant scales. Full simulations of LSST per- 
formance are in progress. To be conservative, we assume 
that (A]) 2 = 10~ 9 , which is roughly the amplitude of a 
z > 2 shear power spectrum at I ~ 1000. Finally, we 
infer from the Sl oan Digital Sky Surv ey galaxy angular 
power spectrum ()Tegmark et al.ll2002f) that the additive 
galaxy power spectrum error due to extinction, photom- 
etry calibration, and seeing will be at (Af ) 2 = 10 -8 level. 

With the assumed systematic uncertainties, we find 
that the joint analysis of LSST BAO and WL can achieve 
~ 1.5% precision on 11 distances from Z2 = 0.29 to z%2 — 
3.6 with weak priors on cosmological parameters. Strong 
priors from Planck will reduce the errors to ~ 0.5% level. 
The corresponding errors on the growth factors are ~ 
3.5% from z\ — 0.14 to z\\ = 3.1 without Planck and 
~ 2% with Planck. These estimates are obtained without 
assuming a specific model for dark energy or modified 
gravity, so they are fairly model independent. 

Finally, we mention that the distance and growth fac- 
tor reconstruction automatically provides a pure met- 
ric measurement of the mean curvature of the universe. 
The joint LSST BAO and WL analysis can achieve 
cr(f2fc) = 0.017 with the above mentioned systematic 
uncertainties and weak priors on the other cosmological 
parameters. Though less impressive than other model- 
dependent forecasts of cr(Ofc) < 10~ 3 (|Knox et alj|2006bt 
IZhan|[2006l ). this estimate does not assume any specific 
dependence of the distance on cosmological parameters. 
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APPENDIX 

INFLUENCE OF COSMOLOGICAL PARAMETERS ON THE DISTANCE AND GROWTH FACTOR CONSTRAINTS 

We check the results with only one cosmological parameter fixed at a time in order to identify which ones contribute 
the most to the difference between the constraints with all the other parameters floating and those with all the 
cosmological parameters fixed in Fig. [T] For BAO distances, the matter density u m , the baryon density tub, and the 
curvature parameter Cl k are the most crucial parameters, and they have comparable impact on the distance errors 
with fl' k weighing more at higher redshift. This is expected qualitatively because uj m and LOb determines the standard 
ruler of BAOs and because the fractional difference between the comoving angular diameter distance Da(z) and the 
comoving distance D(z), 

is larger at higher redshift. The single most important parameter for the BAO growth factors (with galaxy bias 
parameters floating) is the normalization of the matter power spectrum, A|j. The reason is that fixing A 2 ^ removes 
its degeneracy with each Gi for any single galaxy power spectrum. 

The WL results are less obvious. The matter density uj m affects the WL distance constraints the most, then the 
primordial spectral index n s , and then the normalization A|j. Fixing these parameters reduces the uncertainty in the 
amplitude and shape of the matter power spectrum, which is helpful. However, the same argument applies to other 
parameters such as the baryon density LOb and the running of the spectral index a s as well, even though u>t, and a s have 
much less impact on the distance constraints. Yet more puzzling is that the WL growth constraints hardly change 
whether the cosmological parameters (especially ApJ are fixed or not. 
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Fig. 9. — Left panel: Derivatives of the auto shear power spectrum P^ (£) (in the redshift bin 0.7 < z p < 1.05) with respect to a subset 
of parameters as labeled: lnaj m , — lnojj, ±n s , a s . and A. 2 R in dashed lines, and — lnZ54 and U1G4 at 24 = 0.67 in solid lines. Right panel: 
Comparison between (thick lines) and (£) (1-4 < z p < 1.75, thin lines) derivatives with respect to lnoj m , ±n s , — I11.D4, and 

lnG4. To show the difference in the slopes, we normalize the PJ2(£) derivatives to the corresponding P^(t) derivatives at I = 2000 where 
the statistical errors are roughly at the minimum. The derivatives with respect to n s reverse sign near I ~ 100 because the matter power 
spectrum is pivoted at k = 0.05 Mpc -1 = £/Da(z) for n s . The pivot in multipole space moves toward higher £ for shear power spectra in 
(or between) higher redshift bins. 



We do not find an intuitive explanation for the behavior of the WL distance constraints. If the derivatives of the 
shear power spectra with respect to two parameters are nearly degenerate, i.e., roughly proportional to each other, 
then the prior on one parameter will affect the constraint on the other parameter significantly. However, this is not 
the case; as demonstrated in the left panel of Fig. [HI none of the derivatives of the auto shear power spectrum P 2 r 2 7 (£) 



(in the redshift bin 0.7 < z p < 1.05) with respect to lnw,, 



and In A 2 R resemble the derivative with respect to 



\11D4 (z4 = 0.67). Adding to the complexity is that the derivatives vary in both shape and amplitude from one power 
spectrum to another (right panel of Fig. [5]) and that the derivative with respect to one distance parameter differs from 
that with respect to another distance parameter (not shown). This is true for the normalized derivatives din P^* 1 / dq a 
as well. 

The fact that the prior on the power spectrum normalization A 2 R affects the WL distance errors corroborates the 
conclusion in ZK06b that the WL technique measures the distance from the amplitudes of the shear power spectra. 
Closer inspection of equation |T]) shows that uj m also alters the amplitude of the shear power spectrum. The effect is 
twofold: (1) the combined factor ui^/a 2 in the window functions scales (up to a constant factor) the overdensity power 
spectrum A 2 (k) to the potential power spectrum A|(fc) that is directly related to lensing, and (2) an increasing oj m 
moves the broadband turn-over of the transfer function toward small scales, giving a small boost to the shear power 
spectra over the multipole range in our forecasts. This is consistent with dP^ /d\nuj m > 2dP^2 / d\n A|j at I > 40 
and, hence, the prior on w m having a greater impact on the WL distance constraints in Fig. [TJ 

In the right panel of Fig. [5J the derivatives of different shear power spectra with respect to the growth factor 
parameter G4 have almost the same shape. It is likely that this characteristic sets apart the growth parameters 
from the rest, so that the WL growth factor constraints are practically not affected by the priors on the cosmological 
parameters in Fig. [1] To see why these derivatives are nearly proportional to each other, we make the approximations 
A s (k; z) ~ G(z)A s (k) and 

dA 2 (k ; z) . ^(,) A , (fc) ^ D(z 



dG r , 



dG n 



0*3 



D A (z) 



where G(z) is the interpolated growth factor, and 5°{z 



is the Dirac delta function. We have then 



v <xA l3 (z m )r 3 G m A 2 s 



d\nG 



D A {z n 



(Al) 



with Aij(z) — H(z)Da(z)W^ (z)W^ (z) showing that the redshift bins have no influence on the shape of the derivatives 



r /dG n 



dPTr 

become visible in dP?? 



Since equation (|A1[) does not integrate the matter power spectrum over the scales, the BAO features 
/dGi at £ ~ 100 in Fig. 1 

In summary, BAO measurements of the (comoving) distance are most affected by uj m and tub, which determine the 
standard ruler of the sound horizon at the last scattering surface, and by Q' k , which relates the underlying comoving 
distance to the angular diameter distance. WL distance measurements are most affected by us m and A|j, which 
effectively determine the amplitude of the lensing potential power spectrum, and by n s , whereas WL measurements of 
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the growth factor do not require strong priors on cosmological parameters. 
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